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Abstract 

In this paper, we discuss the large-time behavior of solution of a simple kinetic 
model of Boltzmann-Maxwell type, such that the temperature is time decreasing 
and/or time increasing. We show that, under the combined effects of the nonlincarity 
and of the time-monotonicity of the temperature, the kinetic model has non trivial 
quasi-stationary states with power law tails. In order to do this we consider a suitable 
asymptotic limit of the model yielding a Fokkcr-Planck equation for the distribution. 
The same idea is applied to investigate the large-time behavior of an elementary 
kinetic model of economy involving both exchanges between agents and increasing 
and/or decreasing of the mean wealth. In this last case, the large-time behavior of 
the solution shows a Pareto power law tail. Numerical results confirm the previous 
analysis. 

Keywords. Granular gases, overpopulated tails, Boltzmann equation, wealth and 
income distributions, Pareto distribution. 
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1 Introduction 

A well-known phenomenon in the large-time behavior of the Boltzmann equation with 
dissipative interactions is the formation of overpopulated tails |BK00l IEB02al IEB02b] . 
Exact results on the behavior of these tails have been obtained for simplified models, in 
particular for a gas of inelastic Maxwell particles. Our goal here is to show that, at least 
for some simplified kinetic model, the formation of overpopulated tails is not only a be- 
havior typical of systems where there is dissipation of the temperature (cooling) , but more 
generally is a consequence of the fact that the temperature is not conserved. One can 
indeed conjecture that the formation of overpopulated tails in a kinetic model depends 
on the breaking of energy conservation. In kinetic theory of rarefied gases, formation of 
overpopulated tails has been first observed for inelastic Maxwell models [EB02a, EB02bJ. 
Inelastic Maxwell models share with elastic Maxwell molecules the property that the colli- 
sion rate in the Boltzmann equation is independent of the relative velocity of the colliding 
pair. These models are of interest for granular fluids in spatially homogeneous states be- 
cause of the mathematical simplifications resulting from a velocity independent collision 
rate. Among others properties, the inelastic Maxwell models exhibit similarity solutions, 
which represent the intermediate asymptotic of a wide class of initial conditions |BCT03| . 
Recently, the study of a dissipative kinetic model obtained by generalizing the classical 
model known as Kac caricature of a Maxwell gas jPT03| . led to new ideas on the mecha- 
nism of the formation of tails. Indeed, in |PT03j connections between the cooling problem 
for the dissipative model and the classical central limit theorem for stable laws of prob- 
ability theory were found. A second point in favor of our conjecture on tails formation 
comes out from some recent applications to economy of one-dimensional kinetic models 
of Maxwell type |S03|. IP041 ICPT04| . The main physical law here is that a strong econ- 
omy produces growth of the mean wealth (which of course is the opposite phenomenon 
to the dissipation). Nevertheless, the kinetic model led to an immediate explanation of 
the formation of Pareto tails |P897j . Having this in mind, in the next Section we study a 
one-dimensional Boltzmann-like equation which is able to describe both dissipation and 
production of energy. This model has been recently considered in [BBLR03] with the 
aim of recovering exact self-similar solutions. The analysis of [BBLR03] . based on the 
possibility to use Fourier transform techniques to investigate properties of the self-similar 
profiles, shows that in many cases there is evidence of algebraic decay of the velocity 
distributions. On the other hand, except in particular cases, no exact results can be 
achieved. To obtain a almost complete description of the large time behavior of the solu- 
tion, we resort to a different approach. After a brief description of the model, in Section 

2 we introduce a suitable asymptotic analysis, which reduces the Boltzmann equation to 
a Fokker-Planck like equation which has an explicitly computable stationary state with 
power-like tails. In Section 3, we show how similar ideas can be fruitfully applied to 
describe the large-time behavior of some elementary kinetic models of an open economy. 
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Here, the underlying Fokker-Planck equation takes the form of a similar one introduced 
recently in [BM001 ICPT04] . The rest of the paper is devoted to the proof of mathematical 
details. Numerical experiments on the Boltzmann models can be found at the end of the 
paper. 

2 Kinetic models and Fokker-Planck asymptotics 

In this section we will study the large-time behavior of solutions to one-dimensional kinetic 
models of Maxwell-Boltzmann type, where the binary interaction between particles obey 
to the law 

v * = pv + qw, w* = qv + pw; p > q > 0. (1) 

The positive constants p and q represent the interacting parameters, namely the portion 
of the pre-collisional velocities (v, w) which generate the post-collisional ones (v*,w*). As 
it will be clear after Subsection 12.21 the choice p > q is natural in mimicking economic 
interactions, so that we will assume it even in molecular dynamics. As a matter of fact, 
the mixing parameters p and q can be exchanged, which corresponds to the exchange of 
post-collision velocities, without any change in the global collision evolution. 

2.1 Nonconservative kinetic models 

Let f(v, t) denote the distribution of particles with velocity v £ M at time t > 0. The 
kinetic model can be easily derived by standard methods of kinetic theory, considering that 
the change in time of f(v,t) depends on a balance between the gain and loss of particles 
with velocity v due to binary collisions. This leads to the following integro-differential 
equation of Boltzmann type [BBLR03] , 

% = j ]R ( " /(*)/(«)) dw (2) 

where (w*,w*) are the pre-collisional velocities that generate the couple (v,w) after the 
interaction. In ([5)) J = p 2 — q 2 is the Jacobian of the transformation of (v, w) into (v*,w*). 
Note that, since we fixed p > q, the Jacobian J is positive and that the unique situation 
corresponding to J = 1 is obtained taking p = 1 and q = for which the collision operator 
vanishes. 

The kinetic equation ([2]) is the analogous of the Boltzmann equation for Maxwell 
molecules |Bo88l ICIP94j . where the collision frequency is assumed to be constant. Also, 
it presents several similarities with the one-dimensional Kac model |Ka59l IMK66] . It is 
well-known to people working in kinetic theory that this simplification allows for a better 
understanding of the qualitative behavior of the solutions. 

Without loss of generality, we can fix the initial density to satisfy 

/ f (v)dv = l; [ vf (v)dv = [ v 2 f {v)dv = l. (3) 
Jm Jm Jm 
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To avoid the presence of the Jacobian, and to study approximation to the collision operator 
it is extremely convenient to write equation ([2]) in weak form. It corresponds to consider, 
for all smooth functions (f>(v), the equation 

d 



<j){v)f{v, t)dv= / f(v)f(w)(<l>(y*) - <f>{v))dvdw. (4) 

R Jr? 



One can alternatively use the symmetric form 

| / f(v)<Mv)dv = \! f(v)f(w) 
at Jm 2 J M 2 

(5) 

(4>(v*) + 4>{w*) — <p(v) — <p(w))dv dw . 

A remarkable fact is that equations and ([5]) can be studied for all values of the mixing 
parameters p and q, including the case p = q, which could not be considered in equation 

©• 

Choosing 4>(v) = v, (respectively (f>(v) = y2 ) shows that 



m(t) = / vf(v, t) dv = m(0) exp {(p + q — l)t} . (6) 
Jr 

Hence, since the initial density /o satisfies ([3]), m(0) = and m(t) = for all t > 0. 
Consequently, 

E(t) = [ v 2 f(v, t) dv = exp { (p 2 +q 2 -l)t}. (7) 
Jr 

Higher order moments can be evaluated recursively, remarking that the integrals J v n f(v, t) 
obey a closed hierarchy of equations [BKOOj . 

Note that the second moment of the solution is not conserved, unless the collision 
parameters satisfy 

p 2 + q 2 = 1. 

If this is not the case, the energy can grow to infinity or decrease to zero, depending on 
the sign of p 2 + q 2 — 1. In both cases, however, stationary solutions of finite energy do not 
exist, and the large-time behavior of the system can at best be described by self-similar 
solutions. The standard way to look for self-similarity is to scale the solution according 
to the role 

g(v,t) = y/Wfif (yy/E$),t) . (8) 

This scaling implies that f v 2 g(v, t) = 1 for all t > 0. Elementary computations show that 
g = g(v, t) satisfies the equation 



-Q t ~ 2 ( p2 + q2 ~ X ) fa ( vg } = J M (j9( v *)9{™*) ~ g{v)g(w)\ dw. (9) 



In weak form, equation ([9]) reads 



d_ 

If 



1 f d 

4>(v)g(v,t) dv - - (p 2 + q 2 - l) / <t>{v) — (vg) dv 
r 1 Jr av 
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d 



dt 



M 



g(v)g(w)((f>(v*) - (/)(v))dvdw. (10) 

Assuming that 4> vanishes at infinity, we can integrate by parts the second integral on the 
right-hand side of (fTUj) to obtain 

)g(v, t)dv + - (p 2 + q 2 - l) / <j>'{v)vg{v) dv = 
1 JTR 

= [ g{v)g(w){<t>{v*)-cj>{v))dvdw. (11) 
Jm 2 

By the collision rule ([T]), 

v* — v = (p — l)v + qw. 
Let us use a second order Taylor expansion of 4>{v*) around v 

<f>(v*) - <f>(v) = ((p - l)v + qw) <f>'(v) + - {{p - l)v + qw) 2 <P"(v), 
where, for some < 6 < 1 

v = 9v* + (1 - 6)v. 
Inserting this expansion in the collision operator, we obtain the equality 



/ g(v)g(w)(c/)(v*) — 4>(v))dvdw = / g(v)g(w) ((p — l)v + qw) 4>'(v)dv dw+ 
Jm 2 Jm 2 

If g(v)g(w)((p-l)v + qw) 2 0"(v)dvdw + R(p,q), (12) 
2 Jm 2 

where 

R(P,q) = l [ (( P -l)v + qw) 2 (<f ) "(v)-cP"(v))g(v)g(w)dvdw. (13) 
2 Jm 2 

Recalling that g(v,t) satisfies ([3]), we can simplify into fj 12 j) to obtain 

g(v)g(w)((p(v*) — 4>(v))dvdw = (p — 1) / vg(v )4>'(v )dv+ 

JR 2 J Ml 



g(v) ((p - l) 2 v 2 + q 2 ) 4>"{v)dv + R(p, q). (14) 



Substituting (|14p into (jlip . and grouping similar terms, we conclude that g(v,t) satisfies 

d r i 



dt 



(j)(v)g{v,t)dv + -((p-l) 2 + q 2 ) / (f>'(v)vg(v)dv 
m z Jm 



Mi 



g(v) ((p - l) 2 v 2 + q 2 ) cp"{v)dv + fl(p, q). (15) 



1 
2 

Hence, if we set 

r = q 2 t, h(v,r) = g(v,t), (16) 



which implies go(v) = ho(v), h(v,r) satisfies 

b(v)h(v,T)dv + i ( (^-^) + 1 ] J <l>'{v)vh(v)dv 



h (v)[ ( — I' v 2 + l\4>"{v)dv + ^R{p,q). (17) 



M 

Suppose now that the remainder in (|17|) is small for small values of the parameter q. 
Then equation (|17p gives the behavior of g(v,t) for large values of time. Moreover, taking 
P = p(q) such that, for a given constant A 

lim P -^ll = A , (18) 
q ^o q 

equation (I17p is well-approximated by the equation (in weak form) 

^- [ 4>(v)h(v,T)dv + \ (A 2 + l) / cp' (v)vh(v) dv = 

I [ h(v) {X 2 v 2 + 1) <j>"{v)dv. (19) 
Equation (I19p is nothing but the weak form of the Fokker-Planck equation 

which has a unique stationary state of unit mass, given by 

. 3 , 1 

1 \ 



M x {v) = c A ( TT ^ J ) (21) 



where 



CA _ K V , (22) 

v^ r /^l + 2A 2N 



2A 2 

Remark 2.1 T/ie derivation of the Fokker-Planck equation Ii20\) presented in this section 
is largely formal. The main objective here was to show that there are regimes of the mixing 
parameters for which we can expect formation of self-similar solutions to the kinetic model 
with overpopulated tails. We postpone the detailed proof and the mathematical technicalities 
to the second part of the paper. 
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Remark 2.2 The conservative case p + q = 1 can be treated likewise. In this case one 
is forced to choose p = a/1 — q 2 , which gives X = as unique possible value. In the limit 
one then obtains the linear Fokker-Planck equation 

dh 1 (d 2 h d . , A 
Note that in this case the stationary solution M(v) is the Maxwell density 



1 



M(v) = -j=e~ v >\ (24) 



for all q < l/\/2. On the contrary, the non conservative cases are characterized by a X 
different from zero, which produces a stationary state with overpopulated tails. Note that 
from (GUp we have c\ — > l/y/2ir as X — > and thus M(v) = lim^o M\{v). 



Remark 2.3 The possibility to pass to the limit in ( flffj) , with A > 0, is restricted to the 
cases p 2 + q 2 < 1 and p 2 + q 2 > 1, but p > 1. In the case p 2 + q 2 > 1, p < 1, it holds 

(1-p) 2 l-p 
q z l+p 

which forces X towards zero as p — > 1. This case, as the conservative one, gives in the 
limit the linear Fokker-Planck equation. Hence, formation of tails is expected in case of 
dissipation of energy, as well in case of production of energy, but only when the mixing 
parameter p > 1 . 

Remark 2.4 In addition to the conservative case, a second one deserves to be mentioned. 
If p = 1 — q, the kinetic models is nothing but the model for granular dissipative collisions 
introduced and studied in IMY93[ \BKOO[ I BMP 00/ as a one-dimensional caricature of the 
Maxwell-Boltzmann equation IBCG00\. \BC03]j . In this case X = — 1, and the stationary 
state is 

2 / 1 



M 'M = ;lTT^J ' (25) 

This solution solves the kinetic equation U0\) . for any value of the parameter q < 1/2. 

Remark 2.5 The asymptotic procedure considered in this section is the analogue of the 
so-called grazing collision limit of the Boltzmann equation JV98a\ \ V98t$J, which relies in 
concentrating the rate functions on collisions which are grazing, so leaving the collisional 
velocities unchanged. It is well-known that in this (conservative) case, while the Boltzmann 
equation changes into the Landau-Fokker-Planck equation, the stationary distribution re- 
mains of Maxwellian type. 
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2.2 Pareto tails in kinetic models of economy 

In this section we show how to extend the asymptotic analysis of the previous section to 
the case in which the kinetic model describes the time evolution of a density f(v, t), which 
now denotes the distribution of wealth v S iR+ among economic agents at time t > 0. 
The collision ([1]) represents now a trade between individuals. For a deep insight into the 
matter, we address the interested reader to [DYOOl IGGPSoSI I1KK981 IPGS031 IPB03j , and 
to the references therein. With the convention f(v, t) = if v < 0, the kinetic model reads 
[CPT041IP04] . 

1 f(v*)f(w*)-f(v)f(w))dw (26) 



dt J M+ \J 

where (v*,w*) G are the post-trade wealths generated by the couple (v,w) after the 
interaction, along the rule (pQ). As before, the jacobian J = p 2 — q 2 . Since the v- variable 
takes values in -ZR+, the collision rules ([T|) lead to a remarkable difference with respect to 
the case treated in the previous section. The pair {v*,w*) of pre-collision variables that 
generate the pair (v, w) is given by 

pv — qw pw — qv 

V * = —J—> W * = —J— 

While in the former case this pair is always admissible ( v * , w* € M) , in the latter we have 
to discard all pairs of pre-collision variables for which v* < or w* < 0. This shows that, 
for any given v € the product f{v*)f{w*) in (I26P is different from zero only on the 
set B = {(q/p)v < w < (p/q)v}. This implies in other words that, if we fix the wealth 
v £ M + as outcome of a single trade, the other outcome w can only lie on the subset B. 

A great simplification is obtained writing equation (|26|) in weak form, where the pres- 
ence of pre-collision wealths is avoided, 

d_ 

~dl: 



<Kv)f(y,t)dv = / f(v,t)f(w,t)((j)(v*) - <p{v))dvdw. (27) 

K + JFti 



Remark 2.6 The role of the energy is now played by the mean m(t) = J vf(v, t) dv. Note 
however that one can think to equation Ii26\) as the analogous of the isotropic form of a 
hard-sphere Boltzmann equation for a density function f(v',t), v ' € M written with respect 
to energy variable v = (v') 2 /2. In this sense it is again the non conservation of the energy 
that will originate the power law tails. 

To look for self-similarity we scale our solution according to 

g(v,t) =m{t)f(m(t)v,t), (28) 

which implies that f vg(v,t) = 1 for all t > 0. Moreover g = g(v,t) satisfies the equation 

4/ H v )9(v,t)dv - (p + q- 1) / (j){v)^- (vg) dv = 
dt Jm+ Jm+ ov 
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g(v)g(w)((j)(v*) - <t>(v))dvdw. (29) 

Performing the same computations of the previous section, and mutatis mutandis we con- 
clude that g(v,t) satisfies 

37 / <t>{v)g{v , t) dv + q I 4>'{v){v - l)g(v)dv = 
m+ Jm 



dt _ 

\2„,2 i „2„„2 



g(v) ((p - 1) V + q 2 w 2 + 2(p - l)gW) 4>"(y)dv + i2(p, g). (30) 

'M 

The form of the remainder R(p, q) is analogous to that of (|13p . It is clear that the correct 
scaling for small values of the parameter q is now 

r = qt, h(v,r) = g{v,t), (31) 

which implies that h(v, r) satisfies the equation 

— f (p(v)h(v, t) dv + f (f/(v)(v — l)h(v) dv = 
dr J M + Jm 

1 f ^^P- l )\.2,m 



h(v)^ '-v^'i^dv + R^q), (32) 

1 JM 1 

where the remainder Ri is given by 

Rx(p,q) = f [ (qw 2 + 2(p-l)vw)<f>"(v)dv + -R(p,q). 
1 JM+ 1 

Let us consider a parameter p = p{q) such that, for a given constant A > 

lim mzJ? = A. (33) 

q^O q 

Then, equation ([32]) is well-approximated by the equation (in weak form) 

4- I 4>{v)h{v,T)dv + I <f/(y)(v-l)h(v)dv = - [ h(v)v 2 <f>" (v)dv . (34) 
dr Jm Jm 2 Jm 

Equation (|34|) is nothing but the weak form of the Fokker-Planck equation 

which admits a unique stationary state of unit mass, given by the T-distribution [BMOCH 
ICPT04| 



where 

2 

M = 1 + j>1. 

This stationary distribution exhibits a Pareto power law tail for large w's. 

Note that this equation is essentially the same Fokker-Planck equation derived from a 
Lotka-Volterra interaction in [BM001 ISo98l IBMRS02j . 

Remark 2.7 The formal analysis shows that the Fokker-Planck equation follows 
from the kinetic model independently of the sign of the quantity p+q—1, which can produce 
exponential growth of wealth (when positive), or exponential dissipation of wealth (when 
negative). Hence, Pareto tails are produced in both situations, as soon as the compatibility 
condition \33\) holds. As discussed in Remark \2.3\ condition l\33\) is always admissible if 
p + q — 1 < 0, while one has to require p > 1 if p + q — 1 >0. This is quite remarkable 
since it shows that this uneven distribution of money which characterizes most western 
economies may not only be produced as the effect of a growing economy but also under 
critical economical circumstances. 

Remark 2.8 The model studied in jSOS^ corresponds to the choice p = 1 — q+e, with e > 0. 
This interaction implies exponential growth of wealth, and convergence of the solution to 
the Fokker-Planck equation if e = e(q) satisfies 

q^O q 

Since the same limit equation is derived within the choice p = 1 — q — e, we are free to 
choose e negative. The particular choice 

e = -2^ + 2q, 

which implies /i = 3/2 and thus A = 4, leads to the stationary state JS03] 

MM = -^^4#' ( 37 ) 

V27T V b ' Z 



which solves the kinetic equation \29) for all values of the scaling parameter q < 1/4. 



3 The Fourier transform of the kinetic equations 

The formal results of Sections 12.11 and 12.21 suggest that, at least in the limit p — > 1 and 
q —> 0, the large-time behavior of the solution to the kinetic model ([9]) is characterized 
by the presence of overpopulated tails. In what follows, we will justify rigorously this 
behavior, at least for a certain domain of the mixing parameters p and q. We start our 
analysis with a detailed study of the Boltzmann model ([2]). 

The initial value problem for this model can be easily studied using its weak form ([3]). 
Let A^o the space of all probability measures in M. + and by 

M a = ^eM : J \v\ a fi(dv) < +oo,a > oj , (38) 
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the space of all Borel probability measures of finite momentum of order a, equipped with 
the topology of the weak convergence of the measures. 

By a weak solution of the initial value problem for equation ([2]), corresponding to 
the initial probability density fo(w) € J\4 a , a > 2 we shall mean any probability density 
/ € C l {JR,Ai a ) satisfying the weak form of the equation 

^ / <j>{v)f{v,t)dv = f f{v)f(w)(<P{v*)-<p{v))dvdw, (39) 

at JM JM 2 
for t > and all smooth functions </>, and such that for all <p 



t-»o 



M JM 



lim / </>(v)f(v,t)dv = <j){v)fo{v)dv. (40) 



In the rest of this section, we shall study the weak form of equation ([2]), with the normal- 
ization conditions ([3]) . It is equivalent to use the Fourier transform of the equation |Bo88| : 

dTm =Q(f,f)(Z,t), (41) 



Of 

where f(£,t) is the Fourier transform of f(x,t), 



/(£,<)= / e-^f(v,t)dv, 
JM 

and 

Q (1 7) (0 = fipOHtf) - /(0/(o). (42) 

The initial conditions ([3]) turn into 

/(0) = l,/'(0) = 0,f(0) = -l, 
/ G C 2 (M). Hence equation ([41 p can be rewritten as 

^§^ + f(C,t) = f( P 0M)- (43) 

Equation ([4^]) is a special case of equation (4.8) considered by Bobylev and Cercignani in 
[BC03j . Consequently, most of their conclusions applies to the present situation as well. 
The main difference here is that the mixing parameters p and q are allowed to assume 
values bigger than 1. 

We introduce a metric on Ai p by 

ifol^M (44) 

s,gm l£r 

Let us write s = m + a, where m is an integer and < a < 1. In order that d s (F, G) 
be finite, it suffices that F and G have the same moments up to order m. 

The norm (|44p has been introduced in [GTW95] to investigate the trend to equilibrium 
of the solutions to the Boltzmann equation for Maxwell molecules. There, the case s = 
2 + a, a > 0, was considered. Further applications of d s can be found in |CGT99( ICCGOO] 
ITV991IUJT02] . 
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3.1 Uniqueness and asymptotic behavior 



We will now study in details the asymptotic behavior of the scaled function g(v,t). 
As briefly discussed before, a related analysis has been performed in the framework of 
the study of self-similar profiles for the Boltzmann equation for Maxwell molecules in 
[BC031 IBCT03] , Likewise, the role of the Fourier distance in the asymptotic study of 
nonconservative kinetic equations has been evidenced in [PT03j . Consequently, part of 
the results presented here fall into the results of [BC031 IPT03] . and could be skipped. 
Nevertheless, for the sake of completeness, we will discuss the point in an exhaustive way. 

The existence of a solution to equation ([2]) can be seen easily using the same methods 
available for the elastic Kac model. In particular, a solution can be expressed as a Wild 
sum |Bo881 ICGT99] . In order to prove uniqueness, we use the method first introduced 
in [GTW95] . Let fi and /jj be two solutions of the Boltzmann equation (|2|), corresponding 
to initial values f\$ and /2,o satisfying conditions ([3]), and /1, f 2 their Fourier transforms. 
Given any positive constant s, with 2 < s < 3, let us suppose in addition that d s (fi j o, /2,o) 
is bounded. Then, it holds 



(45) 



o [fx - h)_ + _ / 2 (£) _ hipOhiqO - UpOUqO 



at \s\> 

Now, since |/i($)| < 1 (|^(^)| < 1), we obtain 

/i(K)/i(?£)-/2(p£)/ 2 (<7£) 



+ l/ 2 (?£)| 



We set 



h(t,o 



< \fM)\ 



p s < sup 



lei- 



fi(q0-f2(q0 



h-h 



(P s + q s 



(46) 



/i(0-/ 2 (0 



The preceding computation shows that 

<(P S + 

Gronwall's lemma proves at once that 



dh , 



(47) 



< exp{(p s + q s - l)t} 1 1 /i 1 1 00 - 



We have 



Theorem 3.1 Let fi(t) and f2{t) be two solutions of the Boltzmann equation (0), cor- 
responding to initial values f±fi and fop satisfying conditions Then, if for some 
2 < s < 3, d s (fifl, /2,o) is bounded, for all times t > 0, 

d s (h(t),f 2 (t)) <exp{(p s + (? s -l)i}4(/i,o,/2,o). (48) 
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In particular, let /q be a nonnegative density satisfying conditions Then, there exists 

a unique weak solution f{t) of the Boltzmann equation, such that /(O) = fo- In case 
p s + q s — 1 < the distance d s is contracting exponentially in time. 

Let us remark that, given a constant a > 0, 

mn l/i(°fl - £Mj _ fl . , 1in 1/iK) ~ &gfl| _ a . d ,f n uq) 

sup r-j- — a sup — — — a d s {fi, f2). (4yj 

Hence, if g(t) represents the solution f(t) scaled by its energy like in (jSJ), 



?(0 = / 



and from (1491) we obtain the bound 



s/W) r 



d s ( 9l (t),g 2 (t)) = sup = 4(/i(t),/ 2 (t))- (50) 

Using (I48p . we finally conclude that, if <?i(t) and g2{t) are two solutions of the scaled 
Boltzmann equation @, corresponding to initial values f\$ and /2,o satisfying conditions 
([3]), then, if 2 < s < 3, for all times t > 0, 

d s ( gi (t),g 2 (t))<exp[[(p s +q s -l)-^(p 2 + q 2 -l)\ t] d 8 (fi, , / 2 ,o)- (51) 

Let us define, for 5 > 0, 

5 p , f/ (5) = p 2+5 + g 2 + 5 - 1 - It* (p 2 + , 2 - 1) . (52) 

Then, the sign of determines the asymptotic behavior of the distance d s {g\{t),g2{t)). 
In particular, if there exists an interval < 8 < 5 in which S Ptq (5) < 0, we can conclude that 
^2+5(51 (t),92(t)) converges exponentially to zero. Note that, by construction, 5 Pj(? (0) = 0, 
and thus mms{S P) q} < 0. The function ()52[) was first considered by Bobylev and Cercignani 
in |BC03j . The sign of S Ptq , however was studied mainly for p = 1 — q, namely the case 
of the dissipative Boltzmann equation. In Figure [1] a numerical evaluation of the region 
where the minimum of the function is negative for p,q € [0, 2] is reported. 

[Figure 1 about here.] 

Remark 3.2 The behavior of S PA (5) when p 2 + q 2 = 1 is clear. In this case in fact, both 
p and q are less than 1, which implies 

S p , q {S) = p 2+S + q 2+5 - K 0, 

for all 5 > 0. We can draw the same conclusion when p 2 + q 2 > 1, while both p < 1 and 
q<l. 
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Consider now the case p 2 + g 2 > 1, with p > 1. In this case, while — 2 -^- (p 2 + q 2 — l) 
decreases linearly, p 2+s increases exponentially, and the sign of S PtQ (5) becomes positive 
for large values of 5. 

If finally p 2 + q 2 < 1, the sign of S ps {5) for large values of 5 is positive, since, while 



p 2+6 +q 2+& _ 1 > _ 1} 



for 



2 + 5 



5 > 



(p2 +q *- 1 )> 1 

2(p 2 + q 2 ) 

l-(p 2 + q 2 ) ' 



The previous remark indicates that in the general case one can at best hope that S p q {8) 
is negative in an interval (0, 5). To show that this is really the case, one has to investigate 
carefully the behavior of S Ptq (5) in a neighborhood of zero. Since the function S p> „(5) is 
convex for 5 > 0, 

d 2 S p , q {5) 



p 2+5 {\ogpf +q 2+6 {\ogqf > 



2+5 1 



d5 2 



and iS P)(? (0) = 0, in all cases where S Pjq (5) is positive for large values of 5, a sufficient 
condition for S Ptq (5) be negative in some interval < S < 5 is that 



dS p , q (5) 



dd 



<5=0 



p 2 logp + q 2 log q - - (p 2 + q 2 - l) < 0. 



Let us discuss before the case p 2 + q 2 < 1. Given A > 0, we introduce a dependence 
between p and q by setting p = 1 — Xq. Since p > q, this relationship is possible only if 
q < 1/(1 + A). Moreover p 2 + q 2 < 1 requires q < (2A)/(1 + A 2 ). Using this, it is immediate 
to show that there is an interval < q < q in which S' p q (0) < 0. We have 



dS p>q (5) 



dd 



1 



5=0 



G(q) = (1 - Xq) 2 log(l - Xq) +q 2 logq-- ((1 - Xq) + q 2 - l) 



Clearly, G(0) = 0. Moreover 

G'{q) = 2q\ogq - 2A(1 - Ac?) log (1 - Xq), 

and 

G"(q) = 2(1 + log?) + 2A 2 (1 + log(l - Xq)) . 

Now G"(q) < in some interval (0, q\), which implies that G'(q) is decreasing in the same 
interval. But, since G'(0) = 0, G'{q) < in the interval (0, g), where q solves 

2q log q - 2A(1 - Xq) log(l - Xq) = 

Consequently, < at least in the same interval. 

Let us now treat the case p 2 + q 2 > 1, with p > 1. Let us set p = 1 + Ag. We have 



dS P , q (5) 
d5 



1 



5=0 



G( 9 ) = (1 + Ag) 2 log(l + Ag) + g 2 log q - - ((1 + Ag) 2 + g 2 - l) . 
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In this case 

G'(q) = 2q log q + 2A(1 + Xq) log(l + Xq), 

and 

G"{q) = 2(1 + logg) + 2A 2 (1 + log(l + Ag)) . 

As before, G"(q) < in some interval (0, q2), which implies that G'(q) is decreasing in the 
same interval. But, since G'(0) = 0, G'(q) < in the interval (0, g), where q now solves 

2glogg + 2A(l + Ag)log(l + Ag) = 

Consequently, G(q) < at least in the same interval. 
We proved 

Lemma 3.3 Let S Ptq (5), 5 > be the function defined by l\52\) . Given a constant X > 0, if 
p 2 + q 2 < 1, let us define p = 1 - Xq . Then, provided q < min {l/(l + A), (2A)/(1 + A 2 )} 
there exists an interval I- = (0, <5_(g)) such that S p ^ q {5) < for 5 G J_. If p 2 + q 2 > 1, 
and p = 1 + Xq there exists an interval I + = (0, 8+(q)) such that S P:q (5) < for 5 € I + . 
In the remaining cases, namely when p 2 + q 2 = 1 or p 2 + q 2 > 1 but p < 1, S Pt q(5) < for 
all 5>0. 

Lemma 13.31 has important consequences both in the behavior of the solution to the 
Boltzmann equation ([9]), and in the limit procedure introduced in Sections 12.11 and 12.21 
The main consequence of the lemma is contained into the following. 

Theorem 3.4 Let gi(t) and g2(t) be two solutions of the Boltzmann equation cor- 
responding to initial values fi o and satisfying conditions Then, there exists a 
constant 5 > such that, if 2 < s < 2 + 5, for all times t > 0, 

d s ( gi (t),g 2 (t)) < exp{-C s t}d s (f lfi ,h fi ). (53) 

The constant C s = — 5 P)9 (s — 2) is strictly positive, and the distance d s is contracting 
exponentially in time. 



3.2 Convergence to self— similarity 

By means of the estimates of Section 13. H we will now discuss the evolution of moments 
for the solution to equation Q. By construction, the second moment of g(v,t) is constant 
in time, and equal to 1 thanks to the normalization conditions ([3]). We can use the 
computations leading to the Fokker-Planck equation (|23p . choosing 4>{v) = \v\ 2+& , where 
for the moment the positive constant 6 < 1. Suppose that the initial density go(v) = fo(v) 
is such that 

v\ 2+s go(v) dv = m$ < oo. (54) 

Then, since the contribution due to the term -§^(vg(v)) can be evaluated integrating by 
parts, 

f \v\ 2+5 -^-(vg(v))dv = -(2 + 5) f \v\ 2+s g(v,t)dv, 
JR av JR 



Tl 
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we obtain 



d f 'v\ 2+5 g(v,t)dv + (2 + 5) p2 + f 1 / \v\ 2+s g(v,t)dv 



dt 



M z JM 



dv dw (\pv + qw\ 2+s - \v\ 2+s ) g(v)g(w) . (55) 

Let us recover a suitable upper bound for the last integral in (|55|) . Given any two constants 
a, b, and < 5 < 1 the following inequality holds 

(\a\ + \b\) 5 <\a\ s + \b\ s . (56) 

Hence, choosing a = p\v\ and b = q\w\, 

\pv + qw\ 2+S < (pv + qw) 2 (p S \v\ & + q 5 \w\ s " 



Substituting into the right-hand side of (|55|) . recalling that the mean value of g is equal 
to zero, and the second moment of g equal to one, gives 



\pv + qw\ 2+S g(v)g(w) dv dw < 

(pv + qw) 2 (p S \v \ d + q s \w\*\ g(v)g(w) dv dw 

(p 2+5 + q 2+8 ) [ \v\ 2+S g(v)dv + {p 2 q 5 + q 2 p S ) [ 
J JR. J JR, 



\v\ s dv. 

Grouping all these inequalities, and recalling the expression of S Ptq (5) given by ([52 
we obtain the differential inequality 



d_ 

dt 



[ \v\ 2+s g(v,t)dv<S M (6) [ \v\ 2+5 g(v,t)dv + B PiS , (57) 
JR JM 

where, by Holder inequality 



B p ,5<p 2 q d + Q 2 P- (58) 

By Lemma [3.31 for any 5 < 5, S Pt g(6) < 0. In this case, inequality ([57]) gives an upper 
bound for the moment, that reads 

\v\ 2+5 g(v,t)dv <m s + p '* <oo. (59) 

M \dp,q\0)\ 

In the case 5 > 3 we can easily iterate our procedure to obtain that any moment of order 
2 + 5, with 5 < 5 which is bounded initially, remains bounded at any subsequent time. 
The only difference now is that the explicit expression of the bound is more and more 
involved. 

If 5 < 5, we can immediately draw conclusions on the large-time convergence of class 
of probability densities {g(v, t)} t>0 , By virtue of Prokhorov theorem (cfr. [L R79| ) the 
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existence of a uniform bound on moments implies that this class is tight, so that any 
sequence {g{v,t n )} n>0 contains an infinite subsequence which converges weakly to some 
probability measure g^. Thanks to our bound on moments, provided 5 < 6, g^ possesses 
moments of order 2 + 5, for < 5 < 5. 

It is now immediate to show that this limit is unique. To this aim, let us consider two 
initial densities /o,i(i>) and fo,2(v) such that, for some < 6 < 5, 

|2+5 x . r„,\ j„, / , f |„,|2+<5 



v\ 2+ fo,i(v) dv < +oo, / \v\ 2+d f 0t2 (v) dv < +oo. 

Then, by Theorem 13.41 the distance d s (fi(t), f2(t)) between the solutions converges expo- 
nentially to zero with respect to time, as soon as 2 < s < 2 + 5. Let now fo(v) possess 
finite moments of order 2 + 5, with < 5 < 5. Thanks to our previous computations on 
moments, for any fixed time T > 0, the corresponding solution f(v, T) has finite moments 
of order 2 + 5. Choosing fo,i(v) = fo(v), and /o,2(^) = f(v, T) shows that d s (f(t), f(t+T)) 
converges exponentially to zero in time. It turns out that the d s -distance between subse- 
quences converges to zero as soon as S P , q (s — 2) < 0. 

We can now show that the limit function goo(v) is a stationary solution to ©. We 
know that if condition (154j) holds, both the solution g(v,t) to equation ([9]) and goo(v) 
have moments of order 2 + 5, with < 5 < 5 uniformly bounded. Hence, for any t > 0, 
proceeding as in the proof of Theorem 13. 1} we obtain 

ds{Q{g{t),g{t)),Q{9oc,goo)) < (p s + q s + 1) d s {g(t), goo ) . (60) 

This implies the weak* convergence of Q{g{t),g(t)) towards Q(goc, g<x)- In particular, 
due to the equivalence among different metrics which metricize the weak* convergence of 
measures |GTW95l ITV99| , if Cq (M) denotes the set of compactly supported continuously 
differentiable functions, endowed with its natural norm || • ||i, for all cp £ Cq(R), 



<f>(v)Q(g(t),g(t))(v)dv -> / 4>(v)Q( goo , goo )(v) dv. (61) 

M JIR 

On the other hand, for all <f> G Cq(M), integration by parts gives 

d (' 

(v) (vg(v,t)) dv = - v<f>'(v)g(v,t)dv. (62) 

r dv y K 

Since |u(/)'(f)| < |f|||0'||i, and the second moment of g(v,t) is equal to unity, the conver- 
gence of d s (g(t),goo) to zero implies 

v(j)'(v)g{v, t) dv — > / vcj)' \v)goo(v) dv. (63) 
n J ir 



Finally, for all </> G C&{R) it holds 

/ 



m 



>goo(v)) - Q(goo,9oo)(v)^ dv = 0. (64) 



d_ 

\ dv 



This shows that g^ is the unique stationary solution to Q. We have 
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Theorem 3.5 Let 5 > be such that S pq (6) < 0, and let goo(v) be the unique stationary 
solution to equation (0J). Let g(v,t) be the weak solution of the Boltzmann equation (0J), 
corresponding to the initial density /q satisfying 

>\ 2+S fo(v)dv < oo. 



Then, g(v, t) satisfies 

\v\ 2+s g(v, t) dv < cs < oo. 



If < 5 < 1 the constant c$ is given by \59\) . Moreover, g(v, t) converges exponentially 
fast in Fourier metric towards goo{v), and the following bound holds 

d2+s(g(t),9oo) < d 2+< 5(/o,5'oo)exp{-|5p i g((5)|t} (65) 

where S p>q (5) is given by if5£j). 

Depending of the values of the mixing parameters p and q, the stationary solution g^ 
can have overpopulated tails. We can easily check the presence of overpopulated tails by 
looking at the singular part of the Fourier transform [EB02a]. Since the Fourier transform 
of g^ satisfies the equation 

- p2+ i~ l t% + 9(0 = (66) 

we set 

M) = l-\i\ 2 + A\i\ 2+s + ... (67) 

which takes into account the fact that g^ satisfies conditions (J3j) - The leading small £- 
behavior of the singular component will reflect an algebraic tail of the velocity distribution. 
Substitution of expression (|67p into (|66p shows that the coefficient of the power |£| 2+<s is 
AS Pt g(5). Thus, the term ^4|£| 2+<s can appear in the expansion of g(£) as soon as 8 is such 
that S Pt q(5) = 0, 5 > 0. In other words, tails in the stationary distributions are present 
in all cases in which there exists a 5 = 5 > such that S pq (6) = 0. Now the answer is 
contained into Lemma [ 



3.3 The grazing collision asymptotics 

The results of the previous section are at the basis of the rigorous derivation of the Fokker- 
Planck asymptotics formally derived in Section [27T1 and [2T2l Suppose that the initial density 
go(v) = fo(y) satisfies condition (|54p . Using a Taylor expansion, we obtain 

\pv + qw\ 2+S - \v\ 2+S = 
1 



(2 + 6)\v\ d v((p - 1> + qw ) + l(l + S)\v\ S ((p - l)v + qw) 2 , 
where, for some < 9 < 1 



9(pv + qw) + (1 - 9)v. 
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Using this into equality (|55|) . one has 

±[ \v\^ s g{v,t)dv + (2 + s / + f~ 1 I \v\ 2+5 g(v,t)dv 
at Jm a Jm 

(2 + 5) / (\v\ S v((p — l)v + qw)) g(v)g(w) dv dw 
Jm 2 v ' 

+ -(2 + <5)(l + <5) [ \v\ 5 ((p-l)v + qw) 2 g(v)g(w)dvdw. 
2 Jm 2 

Since the momentum of g is equal to zero, we can rewrite (|69p as 



d_ 
dt 



[ \v\ 2 + S g(v,t)dv+ 2 -±^[(p-l) 2 +q 2 ] [ 
Jm 1 Jm 



v\ 2+s g(v,t)dv < 

M 



2 

Assuming < 6 < 1, 

\v\ < {l+p) 5 \v\ 5 + q 5 \w\ S . 
Hence, if \p — l\/q = A, we obtain the bound 



\v\ z+0 g(v,t)dv < 



±[ \v\ 2 + S g(v,t)dv+ 2 -±^q 2 [l + \ 2 } [ 
ai Jm z Jm 

+ i(2 + S)(l + 5)q 2 J 2 ((1 +pf\v\ S + q s \w\ s } (Xv + w) 2 g(v)g{w) dv dw, 



or, what is the same, 

d 
dt 



d_ 

d^ 



(69) 



f-(2 + $)(l + <5) / (^((p-lW + g^V^MdW-u;. (70) 

JJR 2 



v\ 2+5 g(v,t)dv <q 2 C(\,q) \v\ 2+S g(v,t) dv . (71) 
7R 



If we now use (|3"Tj) . it holds 

~ ) h(v,T)dv <C(\,q) [ \v\ 2+s h(v,T)dv , (72) 



namely the uniform boundedness of the (2 + 5)-moment of h(v, r) with respect to q, for 
any fixed time r. 

Consider now the remainder f!13|) . which can be rewritten as 



2 

We need the following 



R(p, q) = — [ (- — -v + w) (<j)"(v) - <j>"(v)) h{v)h{w)dv dw. (73) 



M- 
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Definition 3.6 Let F S (]R), be the class of all real functions <j) on M such that <fi( m \v) is 
Holder continuous of order 5, 

\\4> {m) \\s = sup 19 7 9 { >l < oo, (74) 

V -L W \V-W\ d 

the integer m and the number < 5 < 1 are such that m + 5 = s, and (j)^ denotes the 
m-th derivative of g. 

If G T S (M), with s = 2 + 5, 

U"(v) - 0»| < U"\\ s \v - v\ S < U"\\ 5 \(p - l)v + qw\ 8 . (75) 



In this case, 



R(p,q)<- U"\\s (^-^v + w) h{v)h(w)dv dw < 

2 J m 2 V 1 J 

„2+<5 n 

±—y\\ 5 C 2 (\q) / \v\ 2+S h(v)dv. (76) 
2 Jib? 

Thanks to the uniform bound on (2 + 5)-moment of h(y, r) , it follows that, for any fixed 
time t > 0, 

lim^R(p,q) = (77) 

as soon as G ^(M), with s = 2 + 5. This implies that the limit equation is the 
Fokker-Planck equation (|19p . We proved 

Theorem 3.7 Let the probability density fo G M a > where a = 2 + 5 for some 5 > 0, and 
let the mixing parameters satisfy 

(P-I) 2 x2 



/or some constant A fixed. Then, as q — > 0, for all (ft G J- S (M), with s = 2 + 5 the 
weak solution to the Boltzmann equation |17[ ) for the scaled density h(v,r) = g(v,t), with 
t = q 2 t converges, up to extraction of a subsequence, to a probability density h{w,r). This 
density is a weak solution of the Fokker-Planck equation \19\ ). 



3.4 A comparison of tails 

The result of Section 13.31 establishes a rigorous connection between the collisional kinetic 
equation ((2|) and the Fokker-Planck equation (|19p . The result of Lemma 13.31 coupled 
with the comment of Remark 12.31 then shows that there exists a link between tails of the 
stationary solution of Fokker-Plank and Boltzmann equations. In fact, one can choose 
A 2 > in Theorem [20] if and only if the mixing parameters p and q satisfy the conditions 
of the aforementioned Lemma [3.31 Since the reckoning of the size of the tails is immediate 
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in the Fokker-Planck case, it would be important to know if one can extract from this 
knowledge information about size of the tails of the Boltzmann equation. 

Since the size of tails in the Boltzmann equation is given by the positive root of the 
equation 

s P , q (S) = o, 

where S PtQ is the function (|52p . we will try to extract information by comparing this 
root with the value of the parameter A that characterizes the tails of the Fokker-Planck 
equation. If p > 1, using a Taylor expansion of <S P)(? (<5), with p = 1 + Xq, we obtain 



S p , q (S) 2 + 5 



[X 2 5-l) + 
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2 + 5 



x (1+^,3? 
7 H ^ A —s 



(78) 



where < q < q. This shows that, in the scaling of Theorem 1201 the positive root 5*(q) 
of S Pt q(5) = converges, as q — > to the value 1/A 2 , which characterizes the tails of the 
Fokker-Planck equation. When A > 0, one can easily argue that 5*(q) < 1/A 2 . In this 
case, in fact, 



S p , q (S) _2 + 5 



[(X 2 6-l)+A], 



where A > if q > 0. Hence 



S P!q (l/X 2 ) _2 + 5 



q< 



A>0, 



(79) 



(80) 



that, by virtue of the convexity properties of «S Pi9 (5) implies 5*(q) < 1/A 2 . 

A weaker information can be extracted when p < 1 while p 2 + q 2 < 1. In this case, 
writing p = 1 — Ag, A > 0, we obtain 



5 Pi9 (5) _ 2 + 5 



q 2 2 
where < q < q. Let us set 



q< 



2 + 5 
BX 



(l + 5)£ A 3f^ 



(81) 



(82) 



1 + A 2 ' 

where B < 2. In fact, when p < 1 Lemma 13.31 implies that there is formation of tails only 
when p and q are such that p 2 + q 2 < 1 , which is equivalent to the condition 

2A 



^1 + A 2 " 

Hence, when q satisfies (j82|) . from (j8T|) we obtain the inequality 

,(1 + 5)5 A 4 



(83) 



S p , q (5) > 2 + 5 



[A 2 5-l) -B- 



1 + A 2 



(84) 



Easy computations then show that, if 5 = rX 2 , with < r < 1, the right-hand side of 
is nonnegative as soon as 

3r(l - r)(l + A 2 ) > B(l + rX 2 ). 
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Hence, the biggest value of B for which the right-hand side of (|84p is nonnegative is 
attained when r = 1/2. In this case, B = 3/4, and S*(q) < 2/A 2 . We can collect the 
previous analysis into the following 

Lemma 3.8 Let the mixing parameters satisfy 

(P-1) 2 



q 2 



A 



for some constant A fixed. Then, ifp > 1 the positive root 5*(q) of the equation S pq (8) = 0, 
characterizing the tails of the Boltzmann equation, satisfies the bound 6*(q) < 1/A 2 . If 
p > 1, and at the same time q satisfies the bound \8%fy with B = 3/4, the positive root 
S*(q) of the equation S Pt q(5) = 0, satisfies the bound 5*(q) < 2/A 2 . 



(85) 



We remark here that, in the case p < 1, setting 5 = 1 we obtain an exact formula for 

S P , q {±) 3 
q 2 2 

Choosing A = 1, we get S pq (l) = 0. This case, that corresponds to the conservation of 
momentum in the Boltzmann equation has tails which are invariant with respect to q (see 
Remark 12. 



{X 2 -l)+\q-\x\ 



3.5 Kinetic models of economy 

The analysis of Sections El 13.14 13.21 and 13.31 can be easily extended to equation (|26j) for 
the wealth distribution. We can in fact resort to the methods introduced for the kinetic 
equation on the whole real line simply setting 

F(v, t) = f(v, t)I(v > 0), v € M, (86) 

where 1(A) is the indicator function of the set A. With this notation, equation (j26[) can 
be rewritten as equation (j2J), 



dF(v) 



Of 

Likewise, the weak form Q27|) reads 

d_ 

dt 



f \jF(v*)F(w*) - F(v)F(w)J dw. 



(87) 



4/ (f>(v)F(v,t)dv = I F(v,t)F(w,t)(cj)(v*) - <p(v))dvdw. (J 

M JM 2 



We recall that the role of the energy is now supplied by the mean m(t) = f vF(v,t) dv. 
To look for self-similarity we scale our solution according to 

G(v,t) =m(t)F(m(t)v,t), (89) 
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which implies that J vG(v, t) = 1 for all t > 0. Hence, without loss of generality, if we fix 
the initial density to satisfy 

/ F (v) dv = l; I vF (v) dv = 1 , (90) 
Jm Jm 

the solution G(v,t) satisfies (|9U|) , Then, the same computations of Section show the 
following 

Theorem 3.9 Let f\{t) and f2(t) be two solutions of the Boltzmann equation A26\). cor- 
responding to initial values f\$ and f 2 ^ satisfying conditions A90\) . Then, if for some 
1 < s < 2, d s (fifi, /2,o) is bounded, for all times t > 0, 

d s (h(t),f 2 (t)) <exp{(p s + q s -l)t}d s (f lfi J 2fi ). (91) 

In particular, let /o be a nonnegative density satisfying conditions |5p- Then, there exists 
a unique weak solution f(t) of the Boltzmann equation, such that f(0) = fo- In case 
p s + q s — 1 < the distance d s is contracting exponentially in time. 

Since by ([89]) 

from (|49p we obtain the bound 

d s ( 9l (t),g 2 (t)) = sup 1 ^^"^ 2 ^ 1 = (J-Ydsih^hit)). (92) 

Using (|9"Tj) . we finally conclude that, if g\(t) and g 2 (t) are two solutions of the scaled 
Boltzmann equation (|26p . corresponding to initial values and /2,o satisfying conditions 
([90D, Then, if 1 < s < 2, for all times i > 0, 

d s ( 9l (t),g 2 (t)) < exp {[(p s + g s - 1) - s{p + g - 1)] t} 4(/i,o, / 2 ,o). (93) 

Let us define, for 5 > 0, 

^,g(i) = P 1+S + q 1+S -l-(l + d)(p + q-l). (94) 

Then, the sign of H Pt q now determines the asymptotic behavior of the distance d s (gi(t), g 2 (t)). 
With few differences, the proof leading to Lemma 13.31 can be repeated, obtaining 

Lemma 3.10 Let TZ Ptq (5), 5 > be the function defined by (j94\ )- Given a constant A > 0, 
if P + 1 < 1; let us define p = 1 — Ay^ . Then, provided q < 1/A 2 there exists an interval 
I- = (0,5^(q)) such that TZ p , q {5) < for 6 € If p + q > 1, and p = 1 + X^/q there 
exists a interval 1+ = (0, 5+(q)) such that TZ Ptq (5) < for 6 G J+. In the remaining cases, 
namely when p + q = 1 or p + q > 1 but p < 1, 1Z p ^ q (5) < /or a// 5 > 0. 

The main consequence of Lemma I3.10I is contained into the following. 
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Theorem 3.11 Let gi(t) and g2(t) be two solutions of the Boltzmann equation \26\) . cor- 
responding to initial values fi t o and f 2 fl satisfying conditions \90}) . Then, there exists a 
constant 5 > such that, ifl<s<l + 5, for all times t > 0, 

d s (gi(t),g 2 (t)) < exp{-C s t}d s (f lfi ,f 2)0 ). (95) 

The constant C s = —1Z Ptq (s — 1) is strictly positive, and the distance d s is contracting 
exponentially in time. 

Existence and uniqueness of the stationary solution to equation (|29p follows along the 
same lines of Section 13.21 The main result is now contained into the following. 

Theorem 3.12 Let 5 > be such that 1Z PA (5) < 0, and let goo(v) be the unique stationary 
solution to equation [29]) . Let g(v, t) be the weak solution of the Boltzmann equation [29]) . 
corresponding to the initial density /q satisfying 

v\ 1+S fo(v) dv < oo. 

M+ 

Then, g(v, t) satisfies 

\v\ 1+s g(v, t) dv < c$ < oo, 



J (r 



for some constant cs depending only on p and q. Moreover, g(v,t) converges exponentially 
fast in Fourier metric towards goo{v), and the following bound holds 

d 1+ s{g(t),goo) < d 1+ s{fo, 9oo) exp {-\1Z p , q {6)\t} (96) 

where 1Z PA {5) is given by (9J$ . 

Depending on the values of the mixing parameters p and q, the stationary solution 
can have overpopulated tails. The Fourier transform of g^ satisfies the equation 

- (P + q ~ + G?(0 = G(pt)G(qO. (97) 

We set 

G(0 = l-iC + A\C\ 1+S + ... (98) 



which takes into account the fact that g^ satisfies conditions (|90p . The leading small £- 
behavior of the singular component will reflect an algebraic tail of the velocity distribution. 
Substitution of expression (|98p into (j97|) shows that the coefficient of the power is 
ATZ Pt q(5). Thus, the term j4|£| 1+<5 can appear in the expansion of G(£) as soon as 5 is 
such that H p q (5) = 0,5 > 0. As before, tails in the stationary distributions are present 
in all cases in which there exists a 5 = 5 > such that 7Z p ^ q (5) = 0. Now the answer is 
contained into Lemma I3.101 

Last, one can justify rigorously the passage to the Fokker-Planck equation (|35|) . 
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Theorem 3.13 Let the probability density /o £ M a , where a = 1 + 5 for some 5 > 0, 
and let the mixing parameters satisfy 

Q 

for some A > fixed. Then, as q — > 0, for all G J- S {1R), with s = 1+5 the weak solution to 
the Boltzmann equation KWj) for the scaled density h(v,r) = g(v,t), with r = qt converges, 
up to extraction of a subsequence, to a probability density h(w,r). This density is a weak 
solution of the Fokker- Planck equation [35\) . 

We finally remark that the discussion of Section \3A\ with minor modifications, can be 
adapted to establish connections between the size of the tails of the kinetic and Fokker- 
Planck models. 

4 Numerical examples 

In this paragraph, we shall compare the self-similar stationary results obtained by using 
Monte Carlo simulation of the kinetic model with the stationary state of the Fokker-Planck 
model. The method we adopted is based on Bird's time counter approach at each time 
step followed by a renormalization procedure according to the self-similar scaling used. 
We refer to [PR01] for more details on the use of Monte Carlo method for Boltzmann 
equations. 

We used TV = 5000 particles and perform several iterations until a stationary state is 
reached. The distribution is then averaged over the next 4000 iterations in order to reduce 
statistical fluctuations. Clearly, due to the slow convergence of the Monte Carlo method 
near the tails, some small fluctuations are still present for large velocities. 

Gaussian behavior 

First we consider the case A = for which the steady state of the Fokker-Planck asymptotic 
is the Gaussian (|24j) . We fix p = 1 so that for q < l/v2 we expect Gaussian behavior also 
in the kinetic model. We report the results obtained for q = 0.4 and q = 0.8 in Figure EJ 

[Figure 2 about here.] 

Formation of power laws 

Next we simulate the formation of power laws for positive A. We take p = 1.2 and q = .4 
which correspond to A = 0.5. Keeping the same value of A we then take q = 0.1 and 
p = 1.05. In Figure [2] we plot the results showing convergence towards the Fokker-Planck 
behavior. 

[Figure 3 about here.] 
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A simple growing economy 

We take the case of a growing economy for p = 1 — q + thus corresponding to the limit 
Fokker-Planck steady state (|36|) with A = 2 and \i = 1. As prescribed from our theoretical 
analysis we observe that the equilibrium distribution converges toward the Fokker-Planck 
limit as q goes to 0, with A fixed. The results are reported in Figures HI 

[Figure 4 about here.] 

5 Conclusions 

In this paper we studied the large-time behavior of a simple one-dimensional kinetic 
model of Maxwell type, in two situations, depending wether the velocity variable can take 
values on M or in the former case describing nonconservative models of kinetic the- 
ory of rarefied gases, the latter elementary kinetic models of open economies. In both 
situations it has been shown that the lack of conservation laws leads to situations in 
which the self-similar solution has overpopulated tails. This is particularly important in 
the case of economy, where elementary explanations of the formation of Pareto tails can 
help to handle more complex models of society wealth distribution, where various other 
factors occur. It would be certainly interesting to extend a similar analysis to more real- 
istic situations. Recently, a kinetic model including market returns has been introduced 
[CPT0 4 ■ ! . While for this model the asymptotic convergence to the Fokker-Planck limit 
can be obtained, the property of creation of overpopulated tails has been shown only by 
numerical simulation. In realistic models, in fact, there is a strong correlation among 
densities, due to the constraint of having non-negative wealths after trades, and this ap- 
pears difficult to treat from a mathematical point of view. A further point deserves to be 
mentioned. Recent studies have shown that, while overpopulated tails seem to be generic 
feature of the non-conservative collision mechanism, in the kinetic theory of the Boltz- 
mann equation power-like tails only occur in the borderline case of Maxwell molecules 
interactions |BC03|, IBCT03|, IEB02al IEB02bj , whereas in general collision dissipative pro- 
cesses have stretched exponential tail behaviors [BGP04J. It could be conjectured that the 
corresponding phenomenon in general kinetic models of economy with wealth-depending 
collision frequency manifest a behavior in the form of a lognormal type distribution. 
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perbolic and Kinetic Equations: Asymptotics, Numerics, Applications" HPRN-CT-2002- 
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1 The white domain represents the region where the minimum of the function 
S Pi g is negative for p, q G [0, 2] 

2 Asymptotic behavior for A = of the Fokker-Planck model and the Boltz- 
mann model with p = 1 and q = 0.4, 0.8. Figure on the right is in loglog-scale 

3 Asymptotic behavior for A = 0.5 of the Fokker-Planck model and the Boltz- 
mann model for p = 1.2, q = 0.4 and p = 1.05, q = 0.1. Figure on the right 
is in loglog-scale 

4 Asymptotic behavior for A = 2 of the Fokker-Planck model and the Boltz- 
mann model for p = 1 — q + 2y / g, q = 0.1 and q = 0.01. Figure on the right 
is in loglog-scale 
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Figure 2: Asymptotic behavior for A = of the Fokker-Planck model and the Boltzmann 
model with p = 1 and q = 0.4, 0.8. Figure on the right is in loglog-scale. 



32 



J.=D.50 



X=0.5D 




Figure 3: Asymptotic behavior for A = 0.5 of the Fokker-Planck model and the Boltzmann 
model for p = 1.2, q = 0.4 and p = 1.05, q = 0.1. Figure on the right is in loglog-scale. 
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Figure 4: Asymptotic behavior for A = 2 of the Fokker-Planck model and the Boltzmann 
model for p = 1 — q + q = 0.1 and q = 0.01. Figure on the right is in loglog-scale. 



34 



